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Abstract. Proton acceleration in nearby blazars can be diagnosed measuring their intense TeV 
7 -ray emission. Flux predictions for 1101-1-384 (Mrk421) and 1219-1-285 (ON231), both strong 
EGRET sources (0.1 — 10 GeV), are obtained from model spectra of unsaturated synchrotron pair 
cascades fitted to publicly available multifrequency data. An experimental effort to confirm the 
predicted emission in the range 1-10 TeV would be of great importance for the problems of the 
origin of cosmic rays, the era of galaxy formation and the cosmological distance scale. 


1. Introduction 


1.1. Relativistic jets 

It came as a big surprise to most of the community when numerous powerful extra- 
galactic 7 -ray sources were discovered with the Energetic Gamma Ray Experiment 
Telescope EGRET onboard the Gompton Gamma-Ray Observatory GGRO (Fichtel 
et al. 1994, von Montigny et al. 1994). The sources were soon identified as blazars, 
a subclass of active galactic nuclei (AGN), which are kwown to be the most power¬ 
ful and rapidly variable nonthermal emitters also in other energy bands (Bregman 
1990). The most intriguing fact about the discovery was that the received energy 
flux from many of these enigmatic objects is dominated by y-rays. The precursor 
experiment GOS-B had already indicated the existence of such y-ray dominated 
extragalactic sources, baptized ‘proton quasars’ by Moffat et al. (1983). The TeV 
detection of Mrk421 by the Whipple group makes it probable that the blazar y-ray 
spectra extend to still higher energies. Although it is already difficult to conceive 
that y-ray sources as powerful as Mrk421 with ~ 10'^^ erg s“^ exist in na¬ 
ture, some blazars detected by EGRET had apparent y-ray luminosities of even 
L.y « 10^® erg s-b 

These new challenging observations can be explained assuming that the y-rays 
are emitted in a lighthouse beam directed towards the observer (McBreen 1979). 
The beaming pattern, and thus the blazar phenomenon, arises from the relativistic 
bulk motion of plasma jets ejected from the vicinity of a supermassive black hole 
(Blandford & Konigl 1979). Relativistic radio jets are commonly observed in flat- 
spectrum radio sources and superluminal expansion of the radio structure indicates 
bulk motion with Lorentz factors yj = 10 — 100 (Begelman et al. 1994). Whereas 
most particles in the jet plasma remain in thermal equilibrium, a small fraction of 
them is accelerated to high energies retaining isotropy in the comoving frame by 
pitch-angle scattering. Most probably, particles gain energy by repeatedly cross¬ 
ing collisionless shock fronts in the supersonic jet flow (Drury 1983, Begelman & 
Kirk 1990). A power law energy distribution indicates the stochastic nature of the 
acceleration mechanism. 
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1.2. The radiation process responsible for 7-ray emission 

The high particle energies and the twenty orders of magnitude in frequency of the 
emitted nonthermal radiation are reminiscent of the properties of cosmic rays in 
our own Galaxy where electrons produce synchrotron radio emission and protons 
with Lorentz factors up to 10^^ generate energetic 7 -rays by pion production p + 
p —>■ 'K° + X and their subsequent decay. Proton-matter interactions can not be 
the radiation mechanism responsible for the variable and highly luminous 7 -ray 
emission from blazars, since this would require unrealistically large amounts of 
target matter in the jets. It would also imply a significant enhancement of light- 
elements due to spallation of heavier nuclei which has not been observed (Baldwin 
et al. 1979). However, it was noticed by Blumenthal (1970) and Colgate (1983) 
that the compact radiation fields in extragalactic radio sources represent by far the 
most important target for relativistic protons (cf. also references in Berezinsky et 
al. 1990). At very high energies protons scatter inelastically by photo-production of 
secondaries, i.e. p+j '!T°+p and p+j e'^ +e~ +p. In contrast to the low-density 
interstellar medium where 7 -rays can escape freely, blazars absorb the energetic 7 - 
rays resulting from the decay of secondaries via pair creation 7 -f 7 —> e“'" -I- e“. As 
a consequence, unsaturated pair cascades develope which attenuate the 7 -ray flux 
above the energy e* where = 1 , but which let 7 -rays with energies e < e* 

escape freely (Mannheim et al. 1991). 

There has been some debate whether inverse-Compton scattering or synchrotron 
radiation is the dominant energy loss process of pairs in the blazar jets. In radio 
jets the observed energy density ratio of photons and magnetic field, which governs 
the ratio of the cooling time scales, is a = UjjuB < 1 at a scale of r = I kpc 
(Harris et al. 1994). Assuming a tangled advected magnetic field, which scales as 
B (X r~^, and noting that the photon energy density scales as u-y oc for a 
constant comoving frame luminosity it follows that a < 1 remains invariant. This 
is in accord with X-ray measurements of radio sources (Biermann et al. 1988). 

Theoretically a < 1 is expected, if the target photons are produced locally in the 
jet plasma. The particles with energy density Mrei produce the radiation with energy 
density Uy < Urei- Presumably, particles are accelerated by scattering off magnetic 
field irregularities with energy density usb < ub, hence Uy < Urei < ub- If, on 
the other hand, the relativistic particle pressure would exceed the magnetic field 
pressure, the particles would break out of their magnetic confinement immediately 
turning-off further acceleration. 

Only if the particles responsible for the 7 -ray emission are accelerated very close 
to the AGN, so that the energy density of the external radiation field exceeds that 
instrinsic to the jet, inverse-Compton scattering becomes the dominant cooling 
mechanism. Sikora et al. (1994) have argued that the re-radiation of a central 
thermal UV spectrum by the broad emission line clouds can produce such a strong 
external radiation field, if the 7 -ray production occurs within the central parsec of 
the AGN powering the 7 -ray blazar. Dermer et al. (1992) have argued that direct 
thermal radiation from an accretion disk represents the majority of target photons. 
It was predicted by the latter model that spectral breaks greater than As = 0.5 
should not occur in the 7 -ray regime - in contrast to the observations. Also required 
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is the presence of a source of thermal target radiation dominating the bolometric 
luminosity, although no indication of this big blue bump spectral component is seen 
in the spectra of several blazars with a near-infrared/optical cut-off (Bregman 1990). 
It is a fundamental question, whether or not such a thermal emission component 
is ubiquitous in AGN. Indeed, the emission component can remain undiscovered 
due to Doppler boosting of the nonthermal luminosity by factors > 10'^ or due to 
obscuration by dust. If, however, obscuration is responsible for hiding the thermal 
source, it would also be hidden for the relativistic particles in the jet. Consequently, 
nonthermal optical outbursts should not occur simultanously with y-ray outbursts. 
Another limitation of these models is that the extreme inverse-Compton losses 
postulated imply that an unknown impulsive electron acceleration mechanism must 
be operating. Magnetic reconnection could be such a process. However, the thickness 
of neutral sheets is probably insufficient to obtain the large potential drops required 
for the great values of the electron Lorentz factors. It is also unclear, whether the 
observed power law electron spectra can be produced in this way. 


1.3. Proton-induced cascades 


In the most simple model for blazar emission it is assumed that steady, conical rel¬ 
ativistic jets of Lorentz factor yj with opening angle <I> « yC^ are viewed at a small 
angle 6 (Blandford & Konigl 1979, Mannheim 1993a). Photons emitted isotropically 
in the comoving frame are blueshifted with the Doppler factor 6 = [y(I — /?cos 9)]~^ 
and the apparent luminosity exceeds the comoving frame luminosity by a large 
margin. The outflow is assumed to advect a tangled magnetic held and relativistic 
particles in equipartition with the magnetic held. Most probably shock acceleration 
compensates adiabatic losses during the expansion. Synchrotron-self-absorption fre¬ 
quency, break frequeny (above which the cooling time scale is shorter than the 
dynamical time scale) and emissivity obey a simple scaling with distance r along 
the jet, such that the superposition of spectra from scales between rmin and r^ax 
yields a hat spectrum Si, oc const, up to the break frequency Ub w 10^^ — 10^^ Hz 
above which the spectrum steepens by one power yielding Si, oc In the near- 
infrared/soft X-ray regime the primary electron synchrotron spectrum turns over 
steeply when the cooling time scale becomes shorter than the acceleration time 
scale. For protons, the break energy is reached at 


\mej I-b 2400^®’*’ 


a 


3-10® 



( 1 ) 


Note that a = Uj/ub refers only to the synchrotron photon density in the radio 
through X-ray range, since higher energy photons are unimportant as a target, 
for inverse-Compton scattering because of the Klein-Nishina cross-section and for 
photo-production because of the threshold energy. The proton break energy is much 
greater than the electron break energy, since the proton energy losses are negligible 
compared to the electron energy losses at the same energy. The proton-induced 
luminosity for 7 p < 7p^b has the value 
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where ij = Up/u^ denotes the energy density ratio of relativistic protons and elec¬ 
trons, respectively. It is very likely that protons do not always reach the extremely 
high energy 7 p = 7 p,b- One limitation arises from the condition that the proton 
Larmor radius must remain smaller than the shock radius which is itself less than 
the jet radius (Bell 1978). Hence it follows that for a universal cosmic ray pro¬ 
ton/electron ratio rj it is the proton maximum energy which determines the 7 -ray 
spectrum. 

The brightest 7 -ray emission is expected from the most compact parts of the 
jet where most of the nonthermal target photons with v > Vh are produced, 
i.e. at the jet length r = (jet radius rb,A « rb^)). At smaller scales r < rb 
synchrotron-self-absorption and rapid dynamical evolution darken the jet, whereas 
at larger scales r > rb the radio flux increases and the infrared photon flux de¬ 
creases. Due to the threshold condition for secondary particle production, e.g. 

Dh > 2 • 10 ^^( 7 p/ 10 ^°)“^ Hz for pions, only the infrared and higher frequency 
photons are relevant as a target for the protons. Indeed, the jet would be 7 -ray 
bright also on large scales (steep spectrum radio sources), if the proton maximum 
energy were allowed to reach up to a value 7 p > 7 p,b at any scale in the jet. 
However, this is impossible due to the Larmor radius constraint and the finite ac¬ 
celeration time scale. The greatest values of the proton maximum energy seem to 
be achieved in the hot spots of radio galaxies where 7 p,max ~ 3 • 10^^ « -OTyp^b 
(Biermann & Strittmatter 1987). Notice, however, that blazars with k, lOOLgyn 
require 7 p = 7 p^b ~ 10^° (?? = 100) in the compact parts of the blazar jet. Some of 
these protons escape after isospin flip p -|- 7 ^ n -|- tt’’' as neutrons . The neutrons 
then suffer /3-decay at a distance of = TnC 7 n « 100 kpc from the blazar. In 
this way an extragalactic proton flux without adiabatic losses is injected into the 
intergalactic medium. In the observer’s frame these particles have a Lorentz factor 
7 ~ <57p.b « 10 ^^ 

Tacitly it is assumed that the jets become radiative at a rather large distance 
from their origin. Shocks are generated as the jet propagates through the external 
steep pressure gradient of the elliptical host galaxy. This could happen at a distance 
as far as 1 kpc from the kinematical center (Sanders 1983). 

The spectrum of the proton-induced cascades is complex, producing broken 
power laws and double-humped spectra between X-ray energies and 100 TeV (Mannheim 
et al. 1991). The internal absorption of 7 -rays in ’proton blazars’ (Mannheim 1993a), 
which determines the precise shape of the cascade spectrum, is difficult to assess, 
since the photon density is highly inhomogenous, e.g. varying from the limb to 
the center of the jet and varying with length along the jet. However, the effec¬ 
tive compactness is strongly constrained once the cascade spectrum is known from 
observations between keV and GeV photon energies. 

Since proton-induced cascades inject electromagnetic energy at ultra-high ener¬ 
gies, it is, in fact, the strongest implication of the proton blazar model that TeV 
7 -ray emission can be regarded as a typical property. This is in marked contrast 
to the leptonic models where physical parameters have to be pushed to their lim¬ 
its. The best candidate sources for TeV detection are 7 -ray sources with a hard 
X-ray spectral index = 0.5 — 0.7 indicating a low comoving frame radiation 
compactness. 
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1.4. Merits of 7 -ray measurements 

External absorption of 7 -rays by the infrared/optical radiation field produced inside 
the blazar host galaxy (Mannheim 1993b) and by the cosmic background radiation 
produced by all galaxies (Gould & Schreder 1966, Stecker et al. 1992, MacMinn & 
Primack 1994) must be taken into account. This limits the number of possible TeV 
sources to the few nearest blazars, but it opens up new and challenging possibil¬ 
ities for probing the era of galaxy formation and the cosmic distance scale, since 
'' 77 (^ 7 ) where nir(e-y) is the near-infrared photon density at the pair 

creation threshold Co = 2(meC^)^/e-y cs 0.5/(e-y/TeV) eV. A prerequisite for mea¬ 
suring either Hq or nir via the exponential cut-off produced by external absorption 
is a theoretical prediction of the unabsorbed 7 -ray spectra. To this end I discuss 
below the 7 -ray spectra of 1219-1-285 and Mrk421 in some detail. The strategy is 
to find combined electron synchrotron spectra and cascade spectra fitting to the 
multifrequency data for the same set of parameters. In this way the spectrum is 
highly constrained yielding firm predictions for the TeV range. 


2. TeV emission from Mrk421 and 1219+285 


Internal absorption of 7 -rays inside the blazar jets is crucial to the shape of the 
emerging radiation spectrum. For a steady jet the optical depth in the comoving 
frame is given by 
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where a ~ (1 + ry) ^/ 3 j 7 j$ « (1 + ry) Above the comoving frame energy e* where 
Tyj(€*) = 1 the cascade spectrum steepens by one power, since 


^ 1^,-1 (4) 

''77 

From the proton blazar model fits shown in Figs. 1,2 obtained for the parameters 
in Tab. 1 one can infer that the turnover energies in the observer’s frame (5e* are 
of the order of TeV, viz. 


e*(obs) 
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It is the length scale r = r\y where the jet becomes optically thin for infrared radia¬ 
tion which regulates the internal absorption to values such that TeV is the typical 
turnover energy. Another typical energy can be found as a direct consequence: Pairs 
produced from the optically thick part of the cascade at e* have Lorentz factors 
7 * = e*/2. The synchrotron photons from these pairs have characteristic energies 


e**{ohs) = 
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Below this energy the unsaturated cascade is optically thin, so that the typical 
X-ray spectral index is ax = 0.5 — 0.7 in marked contrast to saturated cascades 
with ax = 0.9 — 1.0 (Svensson 1987). 

No definite predictions for the proton blazar spectra during flux outbursts can be 
made at present, since this would require to solve non-stationary cascade equations 
and to couple the radiative transport with a dynamical model for shock propagation 
and particle acceleration. Most probably, outbursts represent shocks propagating 
into the optically thin part of the jet. The optical outbursts (primary electrons) and 
the hard X-/y-ray outbursts (cascade pairs) should be correlated, since the cascades 
develope in the synchrotron photon target produced by the primary electrons. The 
high-energy cascade emission therefore follows the optical outburst on the very 
short proton photo-production cooling time scale. This time scale is larger than 
the optical outburst decay time scale by a factor of « rj/R where R = L^/Lo 
denotes the 7 /optical luminosity ratio, cf. Eq.(2). If there is a universal value of the 
proton/electron ratio rj, it is predicted that the sources with weakest 7 -ray emission 
compared to the optical emission display the largest delays in 7 -rays. 


TABLE I 

Physical parameters for the proton blazar model fits shown in Figs. 1,2. See Mannheim 
(1993a) for further explanations 



1219-t285 

Mrk421 

B [G] 

4 

40 

V 

30 

100 

rbx [cm] 

7 ■ 10^® 

2 • 10^® 

7 p 

2 ■ 10® 

6-10'^ 

<5 

7 

31 

7j 

5 

20 

e [deg] 

7 

1.5 

L44 

30 

18 

7p/7p,b 

3 • 10“^ 

5•10-^ 


3. Discussion and conclusions 

Abundaning geocentrism it is hard to conveive that baryonic cosmic rays should 
not be present in radio jets. Requiring the Larmor radii of protons to be smaller 
than the jet radius yields a maximum energy of up to 10^^ GeV - large enough 
to make extragalactic radio jets the source of cosmic rays above EeV (Rachen & 
Biermann 1993). 

Recent 7 -ray measurements in the energy band 0.1 — 10 GeV have established 
that the brightest and most compact extragalactic jets are indeed powerful cos¬ 
mic particle accelerators with their maximum luminosity temporarily in the 7 -ray 
band. Proton-induced unsaturated synchrotron cascades can be responsible for the 
observed 7 -ray emission. In marked contrast to the relativistic protons, for which 
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Fig. 1. Multifrequency spectrum of 1219-1-285. Data were obtained from the NED, von 
Montigny et al. (1994), Fink (priv.com.) and the Whipple group (priv.com.). The solid 
line shows the proton blazar model fit for the parameter values in Tab. 1. The dashed line 
shows the expected flux accounting for external absorption (Stecker et al. 1992) 


energy losses become important only near their maximum energy, electrons suffer 
severe energy losses during the entire acceleration process which hinders them to 
reach the high Lorentz factors needed to produce ample TeV photons by inverse- 
Compton scattering of UV-photons. Therefore, TeV 7 -ray emission can be consid¬ 
ered as a typical indicator of proton acceleration in radio jets. 

However, nature has limited our view of the universe by a very restricted horizon 
for 7 -ray eyes. Due to pair-absorption of 7 -rays on cosmic background radiation 
fields the TeV world ends at a luminosity distance of roughly 400 Mpc. Only the few 
nearest blazars are therefore potential TeV sources. Nevertheless, the precise value 
of the TeV horizon depends on two quantities of superior astrophysical interest: The 
Hubble-constant and the near-infrared photon density which entails the conditions 
at birth of normal galaxies. We may hope to shed more light on these cosmological 
cornerstones by measuring and analysing the TeV 7 -ray spectra of nearby blazars. 

Neglecting the non-simultanity of the multifrequency data and on the basic 
assumptions of the proton blazar model (biconical relativistic outflow, isotropy 
of ultra-relativistic particles in the comoving frame, equipartition of relativistic 
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Fig. 2. Multifrequency spectrum of Mrk421 obtained from the NED, von Montigny et al. 
(1994), Punch et al. (1992), Thomas (priv.com.), Kiihn (1994) and Karle (1994). The solid 
line shows the proton blazar model ht for the parameter values in Tab. 1. The dashed line 
shows the expected flux accounting for external absorption (Stecker et al. 1992) 


particle energy and magnetic field energy) I was able to obtain flux predictions in the 
air Cherenkov energy range and above for two nearby blazars Mrk421 (z = 0.308) 
and 1218+285 (z = 0.102). 

In the TeV range internal absorption of y-rays producing a break in the power- 
law spectrum competes with external absorption producing an exponential cut-off. 
The observed TeV y-ray emission from Mrk421 is consistent with little external 
absorption, whereas significant external absorption is likely to be present in the 
spectrum of 1219+285. If the Hubble-constant had a value Ho > 75 km/s/Mpc 
or if the near-infrared photon density had a value less than the expected value 
n.y(0.5eV) « 2 • 10“^ cm“^, the proton-induced cascade flux would have been above 
the Whipple upper limit for this source. However, since external absorption does 
not decrease the observed flux in the air Cherenkov energy range by a large factor, 
1219+285 is very likely to be the next blazar detected by a Cherenkov telescope. 
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Abstract. Proton acceleratiorr irr rrearby blazars carr be diagrrosed measurirrg their irrterrse 
TeV 7 -ray emissiorr. Flux predictions for 1101-|-384 (Mrk421) and 1219-|-285 (ON231), both 
strong EGRET sources (0.1 — 10 GeV), are obtained from model spectra of unsaturated 
synchrotron pair cascades htted to publicly available multifrequency data. An experimen¬ 
tal effort to conhrm the predicted emission in the range 1-10 TeV would be of great 
importance for the problems of the origin of cosmic rays, the era of galaxy formation and 
the cosmological distance scale. 


1. Introduction 


1.1. Relativistic jets 

It came as a big surprise to most of the community when numerous powerful 
extragalactic 7 -ray sources were discovered with the Energetic Gamma Ray 
Experiment Telescope EGRET onboard the Compton Gamma-Ray Obser¬ 
vatory CGRO (Eichtel et ah 1994, von Montigny et ah 1994). The sources 
were soon identihed as blazars, a subclass of active galactic nuclei (AGN), 
which are kwown to be the most powerful and rapidly variable nonthermal 
emitters also in other energy bands (Bregman 1990). The most intriguing 
fact about the discovery was that the received energy flux from many of 
these enigmatic objects is dominated by 7 -rays. The precursor experiment 
COS-B had already indicated the existence of such 7 -ray dominated extra- 
galactic sources, baptized ‘proton quasars’ by Moffat et ah (1983). The TeV 
detection of Mrk421 by the Whipple group makes it probable that the blazar 
7 -ray spectra extend to still higher energies. Although it is already difficult 
to conceive that 7 -ray sources as powerful as Mrk421 with ~ 10“^^ erg 
s“^ exist in nature, some blazars detected by EGRET had apparent 7 -ray 
luminosities of even ~ 10 “^® erg s“^. 

These new challenging observations can be explained assuming that the 7 - 
rays are emitted in a lighthouse beam directed towards the observer (McBreen 
1979). The beaming pattern, and thus the blazar phenomenon, arises from 
the relativistic bulk motion of plasma jets ejected from the vicinity of a 
supermassive black hole (Blandford & Konigl 1979). Relativistic radio jets 
are commonly observed in flat-spectrum radio sources and superluminal 
expansion of the radio structure indicates bulk motion with Lorentz fac¬ 
tors 7 j = 10 — 100 (Begelman et ah 1994). Whereas most particles in the jet 
plasma remain in thermal equilibrium, a small fraction of them is accelerat- 
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ed to high energies retaining isotropy in the comoving frame by pitch-angle 
scattering. Most probably, particles gain energy by repeatedly crossing col¬ 
lisionless shock fronts in the supersonic jet flow (Drury 1983, Begelman & 
Kirk 1990). A power law energy distribution indicates the stochastic nature 
of the acceleration mechanism. 

1.2. The radiation process responsible for 7-ray emission 

The high particle energies and the twenty orders of magnitude in frequency of 
the emitted nonthermal radiation are reminiscent of the properties of cosmic 
rays in our own Galaxy where electrons produce synchrotron radio emission 
and protons with Lorentz factors up to 10^^ generate energetic y-rays by 
pion production p-\-p 7r° -|-A and their subsequent decay. Proton-matter 

interactions can not be the radiation mechanism responsible for the variable 
and highly luminous y-ray emission from blazars, since this would require 
unrealistically large amounts of target matter in the jets. It would also imply 
a signihcant enhancement of light-elements due to spallation of heavier nuclei 
which has not been observed (Baldwin et al. 1979). However, it was noticed 
by Blumenthal (1970) and Colgate (1983) that the compact radiation helds 
in extragalactic radio sources represent by far the most important target for 
relativistic protons (cf. also references in Berezinsky et al. 1990). At very high 
energies protons scatter inelastically by photo-production of secondaries, 
i.e. p -|- y ^ 7r° -|- p and p -|- y ^ e"*" -|- e“ -|- p. In contrast to the low- 
density interstellar medium where y-rays can escape freely, blazars absorb 
the energetic y-rays resulting from the decay of secondaries via pair creation 
y-|-y ^ e“’“-|-e“. As a consequence, unsaturated pair cascades develope which 
attenuate the y-ray flux above the energy e* where = 1, but which 

let y-rays with energies e < e* escape freely (Mannheim et al. 1991). 

There has been some debate whether inverse-Compton scattering or syn¬ 
chrotron radiation is the dominant energy loss process of pairs in the blazar 
jets. In radio jets the observed energy density ratio of photons and magnetic 
held, which governs the ratio of the cooling time scales, is a = u^/ub < 1 
at a scale of r = 1 kpc (Harris et al. 1994). Assuming a tangled advert¬ 
ed magnetic held, which scales as H oc r“^, and noting that the photon 
energy density scales as oc r~^ for a constant comoving frame luminos¬ 
ity it follows that a < 1 remains invariant. This is in accord with X-ray 
measurements of radio sources (Biermann et al. 1988). 

Theoretically a < 1 is expected, if the target photons are produced locally 
in the jet plasma. The particles with energy density Uj-^i produce the radi¬ 
ation with energy density Presumably, particles are accelerated 

by scattering off magnetic held irregularities with energy density ugs < ub , 
hence < Uj-^i < ub- H, on the other hand, the relativistic particle pressure 
would exceed the magnetic held pressure, the particles would break out of 
their magnetic conhnement immediately turning-off further acceleration. 
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Only if the particles responsible for the 7 -ray emission are accelerated 
very close to the AGN, so that the energy density of the external radiation 
held exceeds that instrinsic to the jet, inverse-Compton scattering becomes 
the dominant cooling mechanism. Sikora et ah (1994) have argued that the 
re-radiation of a central thermal UV spectrum by the broad emission line 
clouds can produce such a strong external radiation held, if the 7 -ray pro¬ 
duction occurs within the central parsec of the AGN powering the 7 -ray 
blazar. Dermer et ah (1992) have argued that direct thermal radiation from 
an accretion disk represents the majority of target photons. It was predicted 
by the latter model that spectral breaks greater than As = 0.5 should not 
occur in the 7 -ray regime - in contrast to the observations. Also required is 
the presence of a source of thermal target radiation dominating the bolomet- 
ric luminosity, although no indication of this big blue bump spectral com¬ 
ponent is seen in the spectra of several blazars with a near-infrared/optical 
cut-off (Bregman 1990). It is a fundamental question, whether or not such 
a thermal emission component is ubiquitous in AGN. Indeed, the emission 
component can remain undiscovered due to Doppler boosting of the nonther- 
mal luminosity by factors > 10“^ or due to obscuration by dust. If, however, 
obscuration is responsible for hiding the thermal source, it would also be hid¬ 
den for the relativistic particles in the jet. Consequently, nonthermal optical 
outbursts should not occur simultanously with 7 -ray outbursts. Another 
limitation of these models is that the extreme inverse-Compton losses pos¬ 
tulated imply that an unknown impulsive electron acceleration mechanism 
must be operating. Magnetic reconnection could be such a process. However, 
the thickness of neutral sheets is probably insufficient to obtain the large 
potential drops required for the great values of the electron Lorentz factors. 
It is also unclear, whether the observed power law electron spectra can be 
produced in this way. 

1.3. Proton-induced cascades 

In the most simple model for blazar emission it is assumed that steady, 
conical relativistic jets of Lorentz factor 7 j with opening angle <i> ~ 7 j“^ 
are viewed at a small angle 9 (Blandford & Konigl 1979, Mannheim 1993a). 
Photons emitted isotropically in the comoving frame are blueshifted with the 
Doppler factor b = [ 7(1 — /3cos0)]“^ and the apparent luminosity exceeds 
the comoving frame luminosity by a large margin. The outflow is assumed to 
advect a tangled magnetic held and relativistic particles in equipartition with 
the magnetic held. Most probably shock acceleration compensates adiabatic 
losses during the expansion. Synchrotron-self-absorption frequency, break 
frequeny (above which the cooling time scale is shorter than the dynamical 
time scale) and emissivity obey a simple scaling with distance r along the 
jet, such that the superposition of spectra from scales between Cmin and Umax 
yields a hat spectrum Nj, oc const, up to the break frequency ry, ~ 10 ^^ — 
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10^^ Hz above which the spectrum steepens by one power yielding S,^ oc . 
In the near-infrared/soft X-ray regime the primary electron synchrotron 
spectrum turns over steeply when the cooling time scale becomes shorter 
than the acceleration time scale. For protons, the break energy is reached at 


7p,b 


mp \ ^ 1 -|- a 

mg/ 1-|- 


3 • 10® 


7e,b 

100 


( 1 ) 


Note that a = u^/ub refers only to the synchrotron photon density in the 
radio through X-ray range, since higher energy photons are unimportant 
as a target, for inverse-Compton scattering because of the Klein-Nishina 
cross-section and for photo-production because of the threshold energy. The 
proton break energy is much greater than the electron break energy, since the 
proton energy losses are negligible compared to the electron energy losses at 
the same energy. The proton-induced luminosity for 7 p < 7p^b has the value 


Cease V ( — ) Csyn (2) 

V7p,b; 

where r] = Up/u^ denotes the energy density ratio of relativistic protons and 
electrons, respectively. It is very likely that protons do not always reach the 
extremely high energy 7 p = 7 p,b- One limitation arises from the condition 
that the proton Larmor radius must remain smaller than the shock radius 
which is itself less than the jet radius (Bell 1978). Hence it follows that for 
a universal cosmic ray proton/electron ratio r] it is the proton maximum 
energy which determines the 7 -ray spectrum. 

The brightest 7 -ray emission is expected from the most compact parts 
of the jet where most of the nonthermal target photons with R > ry, are 
produced, i.e. at the jet length r = rn (jet radius rb,_L ~ "rb^h). At smaller 
scales r < Cb synchrotron-self-absorption and rapid dynamical evolution 
darken the jet, whereas at larger scales r > rn the radio flux increases 
and the infrared photon flux decreases. Due to the threshold condition for 
secondary particle production, e.g. //th 7 2 • 10 ^^( 7 p/ 10 ^®)“^ Hz for pions, 
only the infrared and higher frequency photons are relevant as a target 
for the protons. Indeed, the jet would be 7 -ray bright also on large scales 
(steep spectrum radio sources), if the proton maximum energy were allowed 
to reach up to a value 7 p > 7 p^b at any scale in the jet. However, this is 
impossible due to the Larmor radius constraint and the hnite acceleration 
time scale. The greatest values of the proton maximum energy seem to be 
achieved in the hot spots of radio galaxies where 7 p,max ~ 3 • 10 ^^ ~ .OTyp^b 
(Biermann & Strittmatter 1987). Notice, however, that blazars with ~ 
lOOTsyn require 7 p = 7 p^b ~ 10 ^® (7 = 100 ) in the compact parts of the 
blazar jet. Some of these protons escape after isospin flip p -|- 7 ^ ra -|- 7 r+ as 
neutrons . The neutrons then suffer /3-decay at a distance of ~ 


tev.tex - Date: January 16, 1995 Time: 14:36 



TEV GAMMA-RAYS FROM PROTON BLAZARS 


5 


100 kpc from the blazar. In this way an extragalactic proton flux without 
adiabatic losses is injected into the intergalactic medium. In the observer’s 
frame these particles have a Lorentz factor 7 ~ ^ 7 p,b ~ 10^^. 

Tacitly it is assumed that the jets become radiative at a rather large dis¬ 
tance from their origin. Shocks are generated as the jet propagates through 
the external steep pressure gradient of the elliptical host galaxy. This could 
happen at a distance as far as 1 kpc from the kinematical center (Sanders 
1983). 

The spectrum of the proton-induced cascades is complex, producing bro¬ 
ken power laws and double-humped spectra between X-ray energies and 
100 TeV (Mannheim et al. 1991). The internal absorption of 7 -rays in ’pro¬ 
ton blazars’ (Mannheim 1993a), which determines the precise shape of the 
cascade spectrum, is difficult to assess, since the photon density is highly 
inhomogenous, e.g. varying from the limb to the center of the jet and vary¬ 
ing with length along the jet. However, the effective compactness is strongly 
constrained once the cascade spectrum is known from observations between 
keV and GeV photon energies. 

Since proton-induced cascades inject electromagnetic energy at ultra-high 
energies, it is, in fact, the strongest implication of the proton blazar model 
that TeV 7 -ray emission can be regarded as a typical property. This is in 
marked contrast to the leptonic models where physical parameters have to 
be pushed to their limits. The best candidate sources for TeV detection are 
7 -ray sources with a hard X-ray spectral index = 0.5 — 0.7 indicating a 
low comoving frame radiation compactness. 

1.4. Merits of 7-RAY measurements 

External absorption of 7 -rays by the infrared/optical radiation held pro¬ 
duced inside the blazar host galaxy (Mannheim 1993b) and by the cosmic 
background radiation produced by all galaxies (Gould & Schreder 1966, 
Stecker et al. 1992, MacMinn & Primack 1994) must be taken into account. 
This limits the number of possible TeV sources to the few nearest blazars, 
but it opens up new and challenging possibilities for probing the era of galaxy 
formation and the cosmic distance scale, since T^^{e^) oc n\,-{e^)Hf^ where 
^ir(G) near-infrared photon density at the pair creation threshold 

Co = 2(mec‘^ f/e^ ~ 0.5/(e.y/TeV) eV. A prerequisite for measuring either 
Ho or ran via the exponential cut-off produced by external absorption is a 
theoretical prediction of the unabsorbed 7 -ray spectra. To this end I dis¬ 
cuss below the 7 -ray spectra of 1219-|-285 and Mrk421 in some detail. The 
strategy is to hud combined electron synchrotron spectra and cascade spec¬ 
tra htting to the multifrequency data for the same set of parameters. In this 
way the spectrum is highly constrained yielding hrm predictions for the TeV 
range. 
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2. TeV emission from Mrk421 and 1219+285 


Internal absorption of 7 -rays inside the blazar jets is crucial to the shape of 
the emerging radiation spectrum. For a steady jet the optical depth in the 
comoving frame is given by 


’ 2{rri^(? Y j 


^syn 


(e)(7^^(e, e^)rb,_Lde ~ 




(TX 


87rln[ec/eb](iReC^)^ 3 


^b,i(3) 


where a ~ (1 -|- 7 ) ^/ 3 j 7 j<i> ~ (1 -|- 7 ) Above the comoving frame energy e* 
where = 1 the cascade spectrum steepens by one power, since 


A =/c 


1 - exp [-+^] 


('^77 ^ 1 ) 


V 77 


( 4 ) 


From the proton blazar model hts shown in Figs. 1,2 obtained for the param¬ 
eters in Tab. 1 one can infer that the turnover energies in the observer’s frame 
Se* are of the order of TeV, viz. 



ln[ec/eb] 

6 


V 

100 



?’b,l 

3 • 10^® cm 


-1 

TeV (5) 


It is the length scale r = rn where the jet becomes optically thin for infrared 
radiation which regulates the internal absorption to values such that TeV is 
the typical turnover energy. Another typical energy can be found as a direct 
consequence: Pairs produced from the optically thick part of the cascade 
at e* have Lorentz factors 7 * = e*/2. The synchrotron photons from these 
pairs have characteristic energies 


e^yobs) = 3-10-^^SBiXf 



B 

10 G 


e*(obs) 

TeV 


MeV( 6 ) 


Below this energy the unsaturated cascade is optically thin, so that the 
typical X-ray spectral index is ax = 0.5 —0.7 in marked contrast to saturated 
cascades with ax = 0.9 — 1.0 (Svensson 1987). 

No dehnite predictions for the proton blazar spectra during flux outbursts 
can be made at present, since this would require to solve non-stationary 
cascade equations and to couple the radiative transport with a dynamical 
model for shock propagation and particle acceleration. Most probably, out¬ 
bursts represent shocks propagating into the optically thin part of the jet. 
The optical outbursts (primary electrons) and the hard X-/ 7 -ray outbursts 
(cascade pairs) should be correlated, since the cascades develope in the syn¬ 
chrotron photon target produced by the primary electrons. The high-energy 
cascade emission therefore follows the optical outburst on the very short 
proton photo-production cooling time scale. This time scale is larger than 
the optical outburst decay time scale by a factor of ~ rj/R where R = 
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denotes the 7 /optical luminosity ratio, cf. Eq.(2). If there is a universal value 
of the proton/electron ratio 7 , it is predicted that the sources with weakest 
7 -ray emission compared to the optical emission display the largest delays 
in 7 -rays. 


TABLE I 

Physical parameters for the proton 
blazar model hts shown in Figs. 1,2. 
See Mannheim (1993a) for farther 
explanations 



1219-1-285 

Mrk421 

B [G] 

4 

40 

n 

30 

100 

Thi. [cm] 

7 • 10^® 

2 • 10^® 

7p 

2 • 10 '’ 

6 • 10^ 

8 

7 

31 

7j 

5 

20 

0 [deg] 

7 

1.5 

La 

30 

18 

7p/7p.b 

3 • 10“^ 

5-10-" 


3. Discussion and conclusions 

Abundaning geocentrism it is hard to conveive that baryonic cosmic rays 
should not be present in radio jets. Requiring the Larmor radii of protons to 
be smaller than the jet radius yields a maximum energy of up to 10^^ GeV 
- large enough to make extragalactic radio jets the source of cosmic rays 
above EeV (Rachen & Biermann 1993). 

Recent 7 -ray measurements in the energy band 0.1 — 10 GeV have estab¬ 
lished that the brightest and most compact extragalactic jets are indeed 
powerful cosmic particle accelerators with their maximum luminosity tem¬ 
porarily in the 7 -ray band. Proton-induced unsaturated synchrotron cas¬ 
cades can be responsible for the observed 7 -ray emission. In marked contrast 
to the relativistic protons, for which energy losses become important only 
near their maximum energy, electrons suffer severe energy losses during the 
entire acceleration process which hinders them to reach the high Lorentz 
factors needed to produce ample TeV photons by inverse-Compton scatter¬ 
ing of UV-photons. Therefore, TeV 7 -ray emission can be considered as a 
typical indicator of proton acceleration in radio jets. 
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Fig. 1. Multifrequency spectrum of 1219+285. Data were obtained from the NED, von 
Montigny et al. (1994), Fink (priv.com.) and the Whipple group (priv.com.). The solid 
line shows the proton blazar model ht for the parameter values in Tab. 1. The dashed line 
shows the expected flux accounting for external absorption (Stecker et al. 1992) 



However, nature has limited our view of the universe by a very restricted 
horizon for y-ray eyes. Due to pair-absorption of y-rays on cosmic back¬ 
ground radiation helds the TeV world ends at a luminosity distance of 
roughly 400 Mpc. Only the few nearest blazars are therefore potential TeV 
sources. Nevertheless, the precise value of the TeV horizon depends on two 
quantities of superior astrophysical interest: The Hubble-constant and the 
near-infrared photon density which entails the conditions at birth of normal 
galaxies. We may hope to shed more light on these cosmological cornerstones 
by measuring and analysing the TeV y-ray spectra of nearby blazars. 

Neglecting the non-simultanity of the multifrequency data and on the 
basic assumptions of the proton blazar model (biconical relativistic outflow, 
isotropy of ultra-relativistic particles in the comoving frame, equipartition 
of relativistic particle energy and magnetic held energy) I was able to obtain 
hux predictions in the air Cherenkov energy range and above for two nearby 
blazars Mrk421 {z = 0.308) and 1218-^285 {z = 0.102). 
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Fig. 2. Multifrequency spectrum of Mrk421 obtained from the NED, von Montigny et 
al. (1994), Punch et al. (1992), Thomas (priv.com.), Kiihn (1994) and Karle (1994). The 
solid line shows the proton blazar model ht for the parameter values in Tab. 1. The dashed 
line shows the expected flux accounting for external absorption (Stecker et al. 1992) 


Mrk421 (z=0.0308) 



.0 -4.0 -2.0 0.0 2.0 4.0 6.0 8.0 10.0 12.0 14.0 

logJE/eV] 


In the TeV range internal absorption of 7 -rays producing a break in the 
power-law spectrum competes with external absorption producing an expo¬ 
nential cut-off. The observed TeV 7 -ray emission from Mrk421 is consistent 
with little external absorption, whereas signihcant external absorption is 
likely to be present in the spectrum of 1219-|-285. If the Hubble-constant 
had a value Ho > 75 km/s/Mpc or if the near-infrared photon density had 
a value less than the expected value ra.y(0.5eV) ~ 2 • 10“^ cm“®, the proton- 
induced cascade flux would have been above the Whipple upper limit for this 
source. However, since external absorption does not decrease the observed 
flux in the air Cherenkov energy range by a large factor, 1219-|-285 is very 
likely to be the next blazar detected by a Cherenkov telescope. 
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